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ABSTRACT 

We follow the ambipolar-diffusion-driven formation and evolution of a fragment in a magnetically 
supported molecular cloud, until a hydrostatic protostellar core forms at its center. This problem was 
formulated in Paper I. We determine the density, velocity and magnetic field as functions of space 
and time, and the contribution of ambipolar diffusion and Ohniic dissipation to the resolution of the 
magnetic flux problem of star formation. The issue of whether the magnetic field ever decouples 
from the (neutral) matter is also addressed. We also find that the electrons do not decouple from 
the field lines before thermal ionization becomes important and recouples the magnetic field to the 
neutral matter. Ohmic dissipation becomes more effective than ambipolar diffusion as a flux reduction 
mechanism only at the highest densities (a few xlO"'^^ cm~^). In the high-density central parts of the 
core, the magnetic field acquires an almost spatially uniform structure, with a value that, at the 
end of the calculation (n„ « 5 x 10^* cm^'^), is found to be in excellent agreement with meteoritic 
measurements of magnetic fields in the protosolar nebula. Outside the hydrostatic protostellar core, a 
concentration of magnetic flux (a "magnetic wall") forms, which gives rise to a magnetic shock. This 
magnetic shock is the precursor of the repeated shocks found by Tassis & Mouschovias (2005b) which 
cause spasmodic accretion onto the hydrostatic core at later times. 

Subject headings: ISM: clouds - ISM: dust - magnetic fields - MHD - stars: formation - shock waves 



1. INTRODUCTION 

In a companion paper (Tassis & Mouschovias 2006a, 
hereafter Paper I) we formulated the problem of the for- 
mation, by ambipolar diffusion, and evolution of a mag- 
netically supercritical molecular cloud fragment, past 
the phase of ion detachment from the magnetic field 
lines (leaving the electrons as the only attached species). 
The formulation can follow the evolution of the core in 
the opaque regime until the central temperature rises 
to about 1000 K. At higher temperatures, thermal ion- 
ization becomes important and the chemical model no 
longer accurately predicts the degree of ionization in the 
core"^. However, by the time thermal ionization becomes 
important and the increase of the degree of ionization 
reattaches the magnetic field to the matter, the ques- 
tion of whether the magnetic field fully detaches from 
the matter at any point in the evolution of a forming 
protostar will already have been answered. 

During the opaque phase of evolution, several impor- 
tant effects are expected to set in. First, if there is 
ever going to be a complete decoupling of the magnetic 
field from the matter, it must occur after the detach- 
ment of the ions but before thermal ionization becomes 
important. (Complete decoupling refers to conditions 
such that the magnetic field has no significant effect on 
the dynamics of the neutral matter, and the motion of 
the neutral matter no longer affects the magnetic field.) 
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Hence, detachment of the last attached species (the elec- 
trons) from the magnetic field is necessary for magnetic 
decoupling to occur. Second, it is likely that a signif- 
icant part of the magnetic-flux reduction necessary for 
the resolution of the magnetic fiux problem of star forma- 
tion will occur in this phase, due to ambipolar diffusion 
or Ohmic dissipation, or both. (Desch & Mouschovias 
2001 showed that ambipolar diffusion is resurrected dur- 
ing those late stages of contraction of the supercriti- 
cal core.) Third, the formation of the "magnetic wall", 
found bv lTassis &: Mouschovias I ()2005al |bl) to be respon- 
sible for the spasmodic accretion onto the formed proto- 
stellar core, should manifest itself at this advanced stage 
of evolution. 

In this paper, we investigate these effects as well as 
the structure of the evolving core for typical values of 
the free parameters, described below. A detailed param- 
eter study is presented in the following paper (Tassis & 
Mouschovias 2006c, hereafter Paper HI). In § [5] we sum- 
marize the physical processes and assumptions relevant 
to the model cloud, and describe the free parameters of 
the model and their values. In § [3]we describe the results 
for the structure and evolution of the core, the late-time 
behavior of the magnetic field, and we discuss the resolu- 
tion of the magnetic flux problem of star formation. The 
main conclusions are summarized and discussed briefly 
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2. THE PHYSICAL MODEL 



The details of the physical processes included in the 
calculations were presented in Paper I. In summary, 
a model cloud is an axially symmetric, nonrotating, 
thin disk embedded in an external medium of con- 
stant thermal pressure Poxt, with the axis of symme- 
try aligned with the z— axis of a cylindrical-polar co- 
ordinate system r, (j), z. We take advantage of pre- 
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vious calculations d Fiedler fc Mouschovias 1 11992| . 119931 : 
iDesch fc Mouschov ias 2001^7 

which show that a typical 
model cloud relaxes rapidly along magnetic field lines 
and achieves balance between thermal-pressure and grav- 
itational forces in this direction. This balance of forces 
is maintained throughout the evolution. Hence, the di- 
mensionality of the MHD equations can be reduced by 
integrating analytically over z, allowing for a much more 
efficient numerical solution, without loss of any physics 
relevant to the problem. 

The magnetic flux tubes are initially loaded with mass 
by specifying a reference state which has a uniform mag- 
netic field Brof and a column density 



Cn.rcf — 



[1 + {r/l..iYf^ 



(1) 



where CTc.rcf and ^rcf are the centr al column den sity and 
radial lengthscale, respectively (jToomrd |1963( ). The 
cloud boundary is placed initially at r = i?o = 5/rGf, 
where the column density falls to less than 1% of its cen- 
tral value. The reference state is isothermal, at a tem- 
perature of 10 K, corresponding to an isothermal sound 
speed C = 0.188 kms'^. 

The system of equations ( [22a- n] in Paper I) governing 
the evolution of the model cloud contains the following 
(dimensionless) free parameters: Poxt, 'ref, Md,cO- The 
parameter Pext, which is the ratio of the external pres- 
sure and the "gravitational pressure" (= irGa'^ of 
matter in the central flux tube of the reference state, 
enters the calculation through the condition of thermal- 
pressure balance across the cloud surface. The typical 
value used in this run is Pcxt = 0.1. The parameter Ij-d, 
which is the dimensionless lengthscale introduced by the 
profile of the column density in the reference state (eq. 
[1] ) , is a measure of the strength of the initial thermal- 
pressure force relative to the gravitational force inside the 
cloud. The value used in this run is Zrof = 5.5tt > 3^/2, 
which makes the cloud significantly thermally supercrit- 
ical. The parameter //d,cO, which is the initial central 
mass-to-flux ratio in units of the critical value for col- 
lapse, has a value /id.co = 0.25, which makes the model 
cloud magnetically subcritical. 

The above values of the free parameters can be ob- 
tained from the following set of initial values of the 
physical quantities characterizing the reference state: 
rin.co = 2600 cm-3; B.^f = 35.26 ^J.G■, = 0.858 pc; 
Rq = 4.3 pc; Mtot = 99 Mq. Parameters stem- 
ming from the microscopic processes (ionization, chem- 
istry, collisions) are very well constrained observation- 
ally and experimentally, and typical values used in the 
calculations presented here ar e given in Appendix A of 
iTassis fc Mouschovias I ()2005aD . 

An initial equilibrium state is calculated by allowing 
the cloud to relax under the assumption of flux-freezing. 
However, because the cloud is initially magnetically sub- 
critical, only mild readjustments of the column density 
and magnetic field proflles of the reference state take 
place before the equilibrium state is reached. 

The opaque phases of the evolution are treated us- 
ing the adiabatic approximation above a certain density 
riopq, which is an additional free parameter in the model. 
In the case of the typical model cloud, nopq = 10^^ cm^"^. 
Detailed studies of the effect of different values of the free 



parameters on the evolution of the model cloud are pre- 
sented in Paper HI. 

The use of an adiabatic equation of state is a limiting 
case, representing the maximum rate at which the tem- 
perature can increase due to gravitational contraction^ 
- no radiative losses are allowed. The other extreme, in 
which radiative losses are assumed to be efficient enough 
for isothermality to be maintained throughout the calcu- 
lation, is presented in Paper III. The results of an exact 
calculation, in which the evolution of the temperature is 
modeled through detailed radiative transfer calculations, 
will necessarily be bracketed by the results of these two 
extreme cases. 

2.1. Flux- Reduction Timescales 

During the late stages of the contraction of a mag- 
netically supercritical core, the mass-to-flux ratio in the 
central flux tubes changes with time, as a result of the 
combined effects of ambipolar diffusion and Ohniic dissi- 
pation. In order to gauge the relative importance of these 
two ffux-reduction mechanisms and to evaluate their con- 
tribution to the resolution of the magnetic flux prob- 
lem of star formation, it is useful to define characteristic 
timescales for the two processes, that can be calculated 
at any timestep from the values of physical quantities 
and their spatial derivatives. 

In the geometry assumed in the present work, the over- 
all timescale in which the local mass-to-flux ratio evolves, 
, is defined by 



- \d<i>B/dt\ ' 

where d^s /dt = 9$ b / dt+v^-'^^ b is the time derivative 
co moving with the ne utrals. Following the procedure 
in [Mouschovias I ()1991[ ) (eq. [30]), we add and subtract 
Vf ■ V$B and we use the fact that Vf (the effective flux 
advection velocity, including any flux redistribution due 
to Ohmic dissipation, see paper I) defines a frame in 
which <I>B is conserved, to find that 

= {Vn-Vi) ■\7<i>B ^ (Vn-Vf)-^ , (3) 

where the last step assumes axial symmetry and no drift 
velocities between species inside the disk along field lines 
(i.e., the z-direction). Using the generalized Ohm's law 
(eq. [17], Paper I) Vf can be expressed as 



"{i4>+ jo,<l>) 
BzO-p 



(4) 



where Vp is the fiux advection velocity and dp the plasma 
conductivity (see paper I for a detailed discussion of their 
definitions). Hence, equation yields 



1 1 ( 



dr 



1 



$B Ba 



(J0 + JO0)- 



dr 



(5) 



In principle, ambipolar diffusion and Ohmic dissipation may 
also contribute to the heating of the fluid. Then the maximum rate 
at which the temperature could increase is determined by the rate 
of compressional heating (due to gravitational contraction) plus 
the heating rate due to magnetic-diffusion processes. Then we find 
on the basis of the simulation described in this paper that heating 
due to ambipolar diffusion and Ohmic dissipation is five orders of 
magnitude smaller than compressional heating. 
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Fig. 1. — (a) Central number density of neutrals, rin.c, normalized to the central number density of the initial equilibrium state 
("n,cO = 2711 cm"-^), as a function of time. The "star" marks the time of formation of the supercritical core. Evolution of central 
quantities of the typical model cloud as functions of the central neutral density enhancement, J^n,c/'^n,cO, where nn,cO = 2711 cm~'^: (b) 
Abundances of species, (c) Ionization rate at the center of the cloud due to all processes (cosmic rays and radioactive decays), (d) Central 
temperature (in K). (e) z— component of the magnetic field, normalized to its value in the initial equilibrium state (-Bz, cO = 36 /iG). (f) 
Exponent k in the relation Bz,c oc nj; j.. 



The contributions of ambipolar diffusion and Ohmic 
dissipation to the overall flux reduction represented in 
the above equations can be separated on the basis of a 
simple geometrical argument. The term "Ohmic dissi- 
pation" refers to reduction of the magnetic flux due to 
disruption of the current. In the assumed geometry, cur- 
rents can only exist in the 0— direction (j^ and van- 
ish identically). Ambipolar diffusion, on the other hand, 
refers to a redistribution of magnetic flux due to a drift 
between the field lines (and any charged species attached 
to them) and the neutral particles. However, in this par- 
ticular geometry, only drifts in the r— direction result to a 
redistribution of mass among different flux tubes. Drifts 
in the (/)— direction can only shift mass within the same 
flux tube, and do not alter the mass-to-flux ratio of the 
particular flux tube. 

So, the first term on the right-hand side of equation (O 
is the inverse of the timescalc of the mass-to-fiux ratio 
change due to ambipolar diffusion, while the second term 
is the inverse of the timescale of the mass-to-fiux ratio 
change due to Ohmic dissipation. Hence, we may write 



— ^ I'^pr Vnr) 
TAD 1*S 



dr 



and 



so that 
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(6) 
(7) 
(8) 



as expected for processes oc curring in parallel 
(jMouschovias fc Paleologou |[l986[ ). 

3. RESULTS 
3.1. Evolution of Central Quantities 

The model cloud would remain in the initial equilib- 
rium state indefinitely if it were not for ambipolar diffu- 
sion. Figure (IT]i) shows the evolution in time of the cen- 
tral number density n^^c, normalized to its initial value, 
"■n.co — 2711 cm^'^. The evolution of the cloud is followed 
through an enhancement of eleven orders of magnitude in 
central density, after which thermal ionization becomes 
important. Until the formation of a supercritical core, 
denoted in the figure by the position of the "star" (at 
"-n.co — 7.3 X lO'* cm^-^, t = 17 Myr), the evolution is con- 
trolled by ambipolar diffusion and is relatively slow. The 
subsequent evolution is dynamic, although significantly 
slower than free-fall. During this phase, the central den- 
sity increases by more than ten orders of magnitude in 
less than 1.5 Myrs. 

Figure ((TId) exhibits the central abundances of all 
species relative to the neutrals, Xs = ns/n^ (where 
s = e, i, g~, g"*", g°), as functions of the central neutral 
density. At low densities, corresponding to either the 
early stages of the evolution of the central region or to 
the cloud envelope, the grains are primarily negatively 
charged, because of the much more effective electron 
(rather than ion) attachment onto the neutral grains. 
This preferential depletion of electrons onto grains causes 
atomic and molecular ions to be slightly more abundant 
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Fig. 2. — Radial profiles of physical quantities at thirteen different times, {to = 0, = 14.37, t2 = 17.9, = 18.35, t4 = 18.42, 
tg = 18.432, = 18.434, h = 18.435, ts = 18.4354, = 18.4358, tio = 18.4362, in = 18.4364, ti2 = 18.4366). (a) Column density, 
normalized to the central column density of the reference state, ^gf = 5-58 X 10~^ gcm"'^. (b) Temperature in K. (c) Central number 
density of neutrals, nn,c, normalized to the central number density of the initial equilibrium state, as in Fig. JT^). (d) Spatial derivative of 
the density profile s = dlnrtn/dlnr. The initial cloud radius is Ro = 4.23 pc. A "star" on a curve, present only after a supercritical core 
forms, marks the instantaneous radius of the supercritical core. 

than electrons. However, at these low densities, both 
electrons and ions are much more abundant than the 
grain species. The fractional abundances of all grain 
species decrease during the early, ambipolar-diffusion- 
controUed phase of the evolution. This is the result of 
the attachment of charged grains to magnetic field lines, 
which are "left behind" as the neutrals contract through 
field lines under the action of their self-gravity. The same 
tendency is exhibited by the neutral grains because neu- 
tral and charged grains are well coupled through inelastic 
collisions. 



Above densities of ~ 10 cm"'', the grains become the 
main charge carriers. At these densities most grains are 
neutral, while at higher densities positive grains become 
almost as abundant as negative grains. Ions are consid- 
erably more abundant than electrons, due to the faster 
depletion of electrons onto grains. 



The abundances of all charged species decrease rapidly 
above a density « 3 x lO'^^ cm~'^, and then more slowly 
again above a density w 6 x 10^^ cm~^. This behavior is 
the result of changes in the ionization rate, C,, shown in 
Figure ([T}:). At densities up to w 6 x 10^^ cm"^, the ion- 
ization is dominated by the contribution of cosmic rays. 
Its value is constant, C = 5 x 10~^^ s~^ for densities up 
to « 3 x 10^^ cm~'^, when the column density exceeds 
~ 100 g cm~^ and the cosmic rays are appreciably at- 
tenuated. Hence, for higher densities the ionization rate 
declines monotonically, until it is dominated by radioac- 
tive decays (mainly of ""^K) above nn « 6 x 10^^ cm~^, 
and reaches a constant value C = 6.9 x 10"^'^ s"^. Over- 
all, between the initial constant value determined by cos- 
mic rays and the high-density radioactive decay plateau, 
the ionization rate decreases by almost 6 orders of mag- 
nitude. This decrease results in a very low degree of 
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Fig. 3. — Radial profiles of physical quantities at different times, as in Fig. l(2)l- (a) z— component of the magnetic field normalized to 
its value in the reference state, -Brcf = 35.26 fiG. (b) r— component of the magnetic field normalized to -Brof- (c) Spatial derivative of the 
z— component of the magnetic field, b = dlnBz/dlnr. (d) Alfven speed normalized to the sound speed of the reference state, C = 0.188 
kms~^. A "star" on a curve, present only after a supercritical core forms, marks the instantaneous radius of the supercritical core. 



ionization at these high-densities. 

The central temperature as a function of central den- 
sity is shown in Figure fHi)- For this model, the isother- 
mal phase {T — 10 K) lasts up to a density rtopq = 
10^^ cm~'^. For higher densities, we assume an adiabatic 
equation of state with adiabatic index equal to 5/3 un- 
til a temperature of 200 K is reached, which occurs at 
a density 9 x 10^^ cm"-^. At higher temperatures, the 
adiabatic index is changed to 7/5 to account for the ex- 
citation of the rotational degrees of freedom of molecular 
hydrogen. A temperature of 10'^ K is reached at a central 
density of 5 x lO^'' cm~'^, and the calculation is stopped 
at this point, because thermal ionization (which is not 
accounted for) becomes important - see also Footnote 1. 

Figure HJ;) displays the central value of the 
z— component of the magnetic field Bz^c normalized to 
its initial value, i?z,cO — 36 /iG, as a function of central 
density. During the ambipolar-diffusion-controUed evo- 
lution (until the formation of a supercritical core, marked 



by the "star"), B^ ^ increases only by a factor of 1.6. 
After that, a phase of near flux-trapping follows, dur- 
ing which Bz c resembles a power-law B^.c oc nJJ ^ with 
K < 0.5 (see Fig. [If). At densities higher than a few 
X lO^cm"^, Bz c increases less rapidly with density, as 
t he magnetic field starts to decouple from the matter 
(jPesch &: Mouschovias I [200l[) . This gradual process is 
abruptly interrupted because of the change in the equa- 
tion of state. However, the overall decrease in k contin- 
ues, and by the end of the run k has decreased to a mean 
value of 0.1. The oscillations exhibited by k after the 
onset of adiabaticity are the result of the system's at- 
tempt to achieve balance between thermal-pressure and 
gravitational forces in the hydrostatic core. For this rea- 
son, the central density and the central magnetic field 
exhibit small oscillations in time, which in turn manifest 
themselves as the "noise" in k.^ 

® Note that the oscillations in k are physical and are not asso- 
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Fig. 4. — Radial profiles of physical quantities at difli'erent times, as in Fig. ((2J. (a) Ratio of thermal-pressure and gravitational forces, 
(b) Ratio of magnetic and gravitational forces, (c) Ratio of thermal-pressure and magnetic forces. A "star" on a curve, present only after 
a supercritical core forms, marks the instantaneous radius of the supercritical core. 



> - 



0.5 
0.0 

0.5 

1.0 

1.5 

-2.0 
-2.5 



0.5 
0.0 

-0.5 

"1.0 

"1.5 

-2.0 
"2.5 




r / Ro 









5 












"0 


5 




"1 





0) 






> 


"1 


5 




"2 







"2 


5 











5 


















-0 


5 


+ 


"1 







"1 


5 


> 








-2 







"2 


5 







5 








u "0 


5 







> "1 


5 


"2 





"2 


5 



10"''l0"''l0"''l0"^10° 
r / Ro 

t 



"10"°1 0^*10" 
r / Ro 



^10° 




lO^^lO^^lO" 
r / Ro 



^10° 







5 








o 






5 

















5 






-2 





"2 


5 



10"^10"''lO"''lO"^10° 
r / Ro 

tl Q tfi tg t,Q 




Fig. 5. — Radial profiles of the velocities of the (a) neutrals; (b) electrons; (c) ions; (d) negatively-charged grains; (e) positively-charged 
grains; and (f) neutral grains. All velocities are normalized to the sound speed of the reference state [C = 0.188 kms~^). Different curves 
correspond to different times, as in Fig. A "star" on a curve, present only after a supercritical core forms, marks the instantaneous 

radius of the supercritical core. 



3.2. Evolution of Spatial Structure of the Model Cloud 

Figures ([2]) - (fT2|) show radial profiles of physical quan- 
tities at different evolutionary times {tg = 0, ti — 14.37, 
t2 = 17.9, t3 = 18.35, ti = 18.42, ig = 18.432, 
te = 18.434, tj = 18.435, tg = 18.4354, tg = 18.4358, 
tio = 18.4362, til = 18.4364, ti2 = 18.4366 - all times 
given in units of Myr). Each curve corresponds to a time 
such that the central density has increased by a factor 
of ten with respect to the previous curve; the first curve 

ciated with the small-amplitude post-shock oscilla tion s appearing 
in the effective flux advection velocity V[ (see Fig. IIOI I. The latter 
are minimized through the use of artificial viscosity. 



corresponds to the equihbrium state and the last curve 
to the time at which the temperature at the center of the 
model cloud reaches 1000 K. 

3.2.1. Denstties and Temperature 

Figure ([2^) displays the column density of neutral par- 
ticles, in units of the central column density of the ref- 
erence state, as a function of radius at the thirteen dif- 
ferent times given above. In the innermost part of the 
core, which is dominated by thermal-pressure forces that 
smooth out all fluctuations, the column density is almost 
uniform. At late times, the column density outside the 
flat inner region exhibits a profile resembling a broken 
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power-law. It steepens appreciably at small radii (but 
always outside the flat inner region) where a hydrostatic 
protostellar core forms (see § 13.2. 3p . 

Figure ([lb) shows the radial dependence of the cloud 
temperature (in K). The temperature remains constant 
at 10 K for all densities below 10^"'^ cm~^. At higher 
densities, it varies adiabatically with density, with an 
adiabatic index equal to 5/3 for T < 200 K and 7/5 
for higher temperatures. The innermost uniform-density 
part of the core also exhibits spatially uniform tempera- 
ture, since T cx p'^~^- Note that, because the core estab- 
lishes a non-varying asymptotic density profile outside 
the pressure-dominated inner region, the radial location 
of the onset of adiabaticity does not propagate outward 
in time. 

Radial profiles of the number density of neutral par- 
ticles, normalized to the initial central number density 
{nn,cO = 2711 cm~^), are exhibited in Figure ([2j:). The 
logarithmic slope of the density, s — dlnnn/dlnr, is 
shown in Figure ^jX}- There are four prominent fea- 
tures exhibited by s. Inside the supercritical core, s ap- 
proaches but does not become equal to —2 in the region 
10"'' < r/Ro < 10-2 (i.e., lO^^cm <r< 10^^ cm). Its 
mean value in this region is s = —1.8. In addition, there 
are three regions of significant steepening of the density 
profile. The first is located just outside the supercrit- 
ical core (r ~ 0.02i?o)j and is a result of the loss of 
magnetic support for the supercritical core, which con- 
tracts dynamically while the supported envelope cannot 
replenish the inward moving mass fast enough to pre- 
vent the steepening of the profile. The second and the 
third are located at smaller radii. As will be discussed 
below (§ [3T31) . the second dip in s (at r ~ 2 x 10-^i?o) 
corresponds to a balancing of the gravitational force by 
magnetic forces, while the third (at r ~ 2 x 10-^i?o) 
corresponds to a balancing of the gravitational force by 
thermal-pressure forces. Note, however, that s never be- 
comes positive, and therefore the density profile remains 
monotonic throughout the evolution. Figure (|2j;) by it- 
self would give the impression that, outside the uniform- 
density central region, there are three different power- law 
regions: immediately beyond the central uniform den- 
sity region, one might think that s « —3; further out, 
in the region 3 x 10-^i?o ^ ^ 10^^ Rq, it appears that 
s —1.5; and beyond the boundary of the supercritical 
core, s w —2.5. However, Figure ([2]i) reveals that the 
profile of the density is considerably more complicated 
than would be deduced by a visual inspection of Figure 

(Et). 

3.2.2. Magnetic Field and Alfven speed 

The behavior of the magnetic field is demonstrated in 
Figure Q. The z— component of the magnetic field is 
shown in Figure ([3^). At early times, and inside the 
supercritical core, Bz follows a pattern similar to that 
of the column and volume densities, with an innermost 
flat region followed by an almost power-law region. The 
logarithmic slope varies slowly between —0.5 and —1 in 
the region 10""' < r/Ro < 10^^ as is seen in Figure ([St), 
which shows b= din Bz/dlnr. Outside the supercritical 
core the profile flattens significantly and is essentially 
uniform. 

At later times, after the onset of adiabaticity in the 
inner part of the collapsing core, the spatial profile of 



Bz steepens significantly. Outside the location of the 
hydrostatic core, a local maximum appears in Bz at 
r « 10-^i?o, indicating a concentration of magnetic flux 
in this region. This behavior is demonstrated more dra- 
matically in Figure ([St), in which the logarithmic slope 
of Bz acquires large negative values in the region where 
the Bz profile steepens, while it crosses zero and becomes 
positive inward of the location of the local maximum. 

Figure ^jp) shows the spatial dependence of the 
r— component of the magnetic field, generated as a re- 
sult of field-line deformation during the collapse. Inside 
the inner flat-density region, where the magnetic field is 
almost uniform, Br declines rapidly to very small val- 
ues and it goes to zero as r ^ 0. The rapid contrac- 
tion outside the hydrostatic core induces a strong field- 
line deformation, and B^ has a sharp local maximum at 
r « 10-^i?o, where its value becomes comparable to, al- 
though still smaller than, that of Bz; i3r,max/Bz = 0.9. 
In fact, Bz > Br at all radii and all times. 

From Figures ([3^) and ^) it thus becomes clear that 
there are two distinct behaviors exhibited by the mag- 
netic field inside and right outside the hydrostatic core. 
Inside the inner flat-density region, the field configura- 
tion becomes almost spatially uniform, while right out- 
side this region a magnetic flux concentration is ob- 
served, with both the z and r components of the field 
exhibiting local maxima. The value of the z— component 
inside the hydrostatic core is 0.14 G, comparable to 
the protos olar magnetic field as derived from mete- 
oritic data ("Lew ''l988V Stacev. Lovering. fc Parrvlll961t 
Eerndon & Rowe 1974). 

Figure ([3jl) displays the Alfven speed normalized to 
the isothermal sound speed of the initial reference state. 
The local maximum at the radius of the hydrostatic core 
is a result of the local maximum in Bz] however the 
sound speed has also increased at these radii as a result 
of the increased temperature. A most important fea- 
ture in the behavior of va is that inside the supercritical 
core its value becomes comparable to (or smaller than) 
that of the sound speed, which provides an explanation 
for the thermaliz ation of linewidths o bserved in molec- 
ular cloud cores (|Baudrv et al. 1 119811 : iMvers &: Benson I 
Il983f l. In the magnetic star formation theory, the ma- 
terial motions responsible for the observed linewidths 
are attributed to long-wavelength, standing Alfven waves 
(see Mouschovias 1987; Mouschovias & Psaltis 1995; 
Mouschovias, Tassis, & Kunz 2006), with a remarkable 
quantitative agreement between theory and observations. 

3.2.3. Forces 

The radial dependence of the ratio of thermal-pressure 
and gravitational forces is shown in Figure (jl^). The 
gravitational force dominates the thermal force every- 
where (since the model cloud is thermally supercritical) , 
except in the adiabatic core at late times, in which the 
thermal-pressure force exactly balances the gravitational 
force. The radius of the hydrostatic core is initially at 
r = 5 x 10-^i?o and by the end of the run it has shrunk 
to r = 2 X IQ-^Rq. 

Figure shows the radial dependence of the ratio 
of the magnetic and gravitational forces. Even inside the 
dynamically contracting supercritical core, the magnetic 
force is an appreciable fraction of the gravitational force. 
Hence the dynamical contraction of the supercritical core 
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is significantly slower than free-fall. The magnetic force 
becomes much smaher than the gravitational force only 
inside the hydrostatic core, where gravity is balanced by 
the thermal-pressure force. The magnetic force has a 
local maximum with respect to the gravitational force at 
the location of the maximum of the magnetic field seen 
in Figure ([3]). This maximum is located at r > 10~^i?o, 
greater than the radius of the hydrostatic core, defined 
as the radius inside which there is a balance between 
thermal-pressure and gravitational forces (seen in Figure 

Figure ([4)3) exhibits the spatial dependence of the ra- 
tio of thermal-pressure and magnetic forces. The mag- 
netic force is greater than the thermal-pressure force 
(and hence provides the dominant opposition to grav- 
ity) everywhere except in the hydrostatic core, where the 
thermal-pressure force dominates the magnetic force by 
many orders of magnitude. Note, however, that just out- 
side the hydrostatic core a local minimum in the thcrmal- 
to-magnetic force ratio forms, caused by the increased 



magnetic field strength shown in Figures ([3^) and ^^). 

3.2.4. Velocities 

The velocity of the neutrals, normalized to the isother- 
mal sound speed of the initial state, is plotted as a func- 
tion of radius in Figure ([5^). Before the onset of adia- 
baticity, the contraction speed is equal to a few tenths of 
the isothermal sound speed at the radius of the supercrit- 
ical core and, at smaller radii, it increases to 0.5 — 1.0 C. 
At yet smaller radii it decreases again and vanishes at 
r — 0. After the innermost parts become adiabatic, 
the mass accumulating in the hydrostatic core causes 
an accelerated infall, and the magnitude of the veloc- 
ity maximum becomes almost equal to 2 C. The in- 
fall speed Vn decreases precipitously inward of the max- 
imum (at r/RQ = 7 X 10~^) and then again inward of 
r = 9 X 10~®i?o, at which a second, local maximum ex- 
ists. This is the result of two shocks at the two locations, 
the outer one being magnetic in origin - it coincides with 
the local maximum of the magnetic force - and the in- 
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ner one being an accretion shock at the boundary of the 
hydrostatic core. 

Comparing the behavior of the velocities of the other 
species with that of the neutrals, we observe that, at the 
late times of interest, the velocity profiles of the ions and 
all grain species (Figs. [SJ; - [Sf ) are almost identical with 
those of the neutrals. Both the thermal and the magnetic 
shocks are experienced by the ions and the grains. The 
situation is different for the electrons (Fig. [SJd), which 
experience only the magnetic shock, while their velocity 
profile does not exhibit any significant feature at the lo- 
cation of the thermal shock. This is because the electrons 
are still well attached to the field lines (see discussion of 
Fig. EJd below) and are decelerated to almost zero veloc- 
ity behind the magnetic shock. The margin for further 
electron deceleration due to the thermal shock (which 
occurs at smaller radii) is therefore much smaller than 
that for the other species. Hence, the electrons appear 
to experience only the outermost, magnetic shock. 

3.2.5. Attachment Parameters 

Figures ([6]) and (O show radial profiles of the attach- 
ment parameters for all species other than the neutrals at 
the thirteen different times. The attachment parameter 
for any species s is defined as 



A, = 



Vf - Vn 



(9) 



where Vf is the effective flux advection velocity (the veloc- 
ity at which the field lines appear to be moving inward). 
If As « I, then species s is attached to the field lines. 
If As ^ 1, then species s is detached from the magnetic 
field lines and its motion follows that of the neutrals. 

At low densities, all charged species and the neutral 
grains (which are well coupled to the charged grains 
through inelastic collisions) are attached to the field 
lines. The grain species begin to detach first (at a radius 
r/Ra w 10^2 ^ X.3 x IO"cm = 8714 AU, and a den- 
sity of rin^c ~ lO^cm"^), just as the magnetically (and 
thermally) supercritical core forms, and have completely 



detached by n^^c ~ 10^ cm The ions begin to detach 
later, at a radius r/Ro « 10"'' = 1.3 x 10^^ cm « 87 AU 
and density nn,c ~ 10^ cm"'^, and have detached com- 
pletely by fia^c ~ 10^°cm~'^. The ions are therefore still 
attached to the field lines at the location of the magnetic 
shock. 

The electrons, as seen in Figure ^}p), have only par- 
tially detached from the field lines by the end of the run, 
when the temperature reaches 1000 K at the center of 
the cloud and thermal ionization is about to become im- 
portant, increase the degree of ionization, and reattach 
the charged species to the field lines. Even at late times 
and in the innermost parts of the hydrostatic core, the 
electron attachment parameter does not become smaller 
than 0.9. A similar conclusion can also be reached by 
examining the direct attachment parameter of the elec- 
trons, the value of which at the center of the cloud is 
plotted in Figure ([7^) as a function of the (normalized) 
central neutral density. It is the case that UeTen ^ 1 
essentially throughout the evolution. 

3.2.6. Contribution to Electric Current Density by Different 

Species 

Figure ^ shows the contribution of different charged 
species to the overall current, as a function of radius at 
different times. At early times, electrons and ions are 
the main current carriers, and their contribution to the 
overall current density is comparable and approximately 
three orders of magnitude greater than that of positive 
and negative grains. At densities high enough so that 
the electron depletion onto the neutral grains causes the 
ion abundance to significantly exceed the abundance of 
the electrons, the ion current becomes greater than the 
electron current by almost a factor of ten. Finally, in 
the innermost region, the current density decays sharply 
due to the increasing effect of Ohmic dissipation (see 
Fig. If 2|1 . In this regime, grains become the main current 
carriers, their contribution being almost two orders of 
magnitude greater than that of the electrons and ions. 
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3.2.7. Mass-to-Flux Ratio and Drift Velocities 

The evolution of the mass-to-flux ratio in the model 
cloud is demonstrated in Figure ([9]). Figure ([9li) shows 
radial profiles of the mass-to-flux ratio, normalized to the 
critical value for collapse, at the thirteen times given in 
§ 13.21 The "star" marks the radius of the supercritical 
core, when the normalized mass-to-flux ratio becomes 
equal to unity. The mass-to-flux ratio is a three-slope 
power law with an inner flat region, corresponding to 
the flat-density, uniform-magnetic-field core. Inside the 
supercritical core but outside the innermost flat region 
it exhibits two distinct behaviors: at outer radii it in- 
creases only mildly with decreasing radius, while at in- 
ner radii, inside the hydrostatic core, it increases sharply, 
and reaches a value a few hundred times greater than the 
critical value by the end of the run. Outside the super- 
critical core, the mass-to-flux ratio decreases with radius 
and is independent of time. 

A similar behavior can be seen in Figure (Od), where 



the normalized mass-to-flux ratio at the center of the 
cloud is plotted as a function of central density en- 
hancement. Initially, during the ambipolar-diffusion- 
controUed phase, the central mass-to-flux ratio increases 
sharply, until it reaches the critical value, and a supercrit- 
ical core forms. During the dynamical contraction stage, 
and while the ions are still attached to the magnetic field 
lines, the mass-to-flux ratio increases only mildly. Fi- 
nally, after the ion detachment, a much steeper increase, 
by approximately two orders of magnitude, takes place. 
This is an important contribution toward the resolution 
of the magnetic flux problem of star formation, as the 
magnetic flux remaining at this stage approaches that ob- 
served in strongly magnetic stars. However, additional 
magnetic flux loss is required to account for the entire 
observed range of magnetic fields in newborn stars. 

As discussed in detail in Paper I, the motion of the field 
lines can be described by an effective flux-advection ve- 
locity, V{, which is the apparent velocity in the lab frame 
with which the field lines are transported inward. When 
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one or more charged species are attached to the field 
Hnes, their velocity is « Vf. When the field completely 
detaches from the matter, — > and the field lines ap- 
pear to be "held in place" as matter drifts inward through 
them. 

The magnetic field lines threading the supercritical core 
extend into the "external" medium (the envelope) and 
the magnetic field is maintained by currents in the par- 
ent cloud envelope. This is well justified because ob- 
servations of the morphology of magnetic fields in and 
around molecular cloud "clumps" show fields of length- 
scale greater than that of the clumps, threading the 
clumps and extending s r noothly into the pareii t cloud 
(|Hildebrand et al. Ifl999l: iMcGreeor e^"anfl99l . Flux 
redistribution precesses (e.g., ambipolar diffusion) and 
dissipative processes (e.g., Ohmic dissipation) cannot 
therefore change the topology of any particular fiux tube, 
or decrease the magnetic-field strength below that im- 
plied by the "external" currents. They only limit the 
speed at which field lines are transported inward, to a 
value smaller than that of the velocity of the neutrals, 
with the consequence that mass can accumulate in the 
central flux tubes, increasing the mass-to-flux ratio there. 
Thus the total magnetic flux of the model cloud (but not 
that of the protostellar cores forming in it) is conserved. 

Radial profiles of Uf are shown in Figure (|10b .V Since, 
as we discussed above, the electrons do not detach from 
the magnetic field lines during the evolution, the radial 
and temporal behavior of Uf resembles that of Vc (see Fig. 
[S}d) , with the prominent deceleration feature correspond- 
ing to the magnetic shock. 

Figure ((TUJs) shows the drift velocity between the field 
lines and the neutrals, Wf — Wn- It exhibits a local maxi- 
mum in the model cloud envelope, where ambipolar dif- 
fusion continues to control the infall of the neutrals, with 
\vi\ <C |wn| and Wn ~ — wd- It then becomes small for a 
significant part of the dynamically contracting supercrit- 
ical core, where the field is almost (although not exactly) 
frozen in the matter. Finally, it becomes significant again 
in the inner parts of the core, where the degree of ioniza- 



tion is very low and the coUisional drag on the neutrals 
decreases. 

3.2.8. Mass and Magnetic Flux 

Figures pTk ) and pTb ) show, respectively, radial pro- 
files of the cumulative mass and magnetic flux at the 
thirteen different times given above. The "star" on a 
curve marks the location of the radius of the supercriti- 
cal core, which contains only a few percent of the cloud 
mass. The mass included inside the hydrostatic core by 
the end of the run is w 10^^ M©. 

At late times, both quantities exhibit a narrow plateau. 
The change from an isothermal to an adiabatic equation 
of state at > 10^^ cm~^ results in the formation of 
a high-density hydrostatic core, corresponding to a cen- 
tral mass concentration at the innermost part of the the 
supercritical core. The plateau in the mass proflle ap- 
pears because outside the hydrostatic core r « 10^^i?o 
the density of the disk declines rapidly (see Fig. [2]) 
and as a result the mass content of the plateau region 
(10^^i?o ^ r < 2 X 10~^i?o) is a very small fraction of 
the mass of the hydrostatic core. 

The deceleration of the infall of field lines at r « 
10~^i?o (Fig. [TUk ) results in a concentration of flux 
outside the hydrostatic core. This concentration dom- 
inates the magnetic flux profile in the plateau region 

3.2.9. Ambipolar Diffusion and Ohmic Dissipation 

The relative effectiveness of ambipolar diffusion and 
Ohmic dissipation in increasing the mass-to-flux ra- 
tio is demonstrated by Figure (HH), which shows the 
ratio of the ambipolar-diffusion and Ohmic-dissipation 
timescales for the central cell as a function of density 
(Fig. \T^). and radial profiles of the same quantity at the 
thirteen different times (Fig. [T^). At low densities and 
large radii, ambipolar diffusion dominates Ohmic dissipa- 
tion. The two processes become equally important only 
at a density 5 x 10^^ cm~'^, while at higher densities 
Ohmic dissipation operates on a timescale shorter than 
that of ambipolar diffusion. 
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Fig. 10. — Radial profile of (a) the effective flux-advection velocity, and (b) the drift velocity (difference between the effective flux- 
advection velocity and the velocity of the neutrals), in units of the sound speed of the reference state. Different curves correspond to 
different times, as in Fig. l[2]l. A "star" on a curve, present only after a supercritical core forms, marks the instantaneous radius of the 
supercritical core. 

1000 K. Still, the electrons are very few and they are 
very ineffective in transmitting the magnetic force to the 
neutrals through collisions. Therefore, the magnetic field 
in the innermost part of the contracting core is dynami- 
cally insignificant by the end of the run. 

The grains begin to detach from the field lines at a 
central neutral density nn,c ~ 10'* cm~'^, at radius « 10'* 
AU. Their detachment is complete by the time the central 
neutral density increases to « 10* cm~'^. By contrast, 
the ions begin to detach when rin^c ~ 10* cm~'^, in the 
region r < 10^ AU, and their detachment is complete 
by the time nn,c ~ 10** cm""^. The electrons in this 
representative model cloud, do not detach from the field 
lines by the time thermal ionization sets in. 



Interestingly, when the grains become the primary 
charge carriers, at densities « 2 x 10*° cm~^, ambipo- 
lar diffusion still dominates Ohmic dissipation as a flux- 
reduction m echan ism, contrary to the expectations of 
iNakano "fc~ Umcbav ashi I (|1986al lbl'). Instead, it is when 
the grains become the main current carriers, at densities 



w 10 cm~'^ that Ohmic dissipation dominates. The 
region in which this occurs coincides with the region in 
which the electric currents diminish (r < 10~^i?o; see 
Fig. El). 

At the location of the magnetic shock, ambipo- 
lar diffusion is always more effective than Ohmic dis- 
sipation. This is c o nsistent with the findings in 
iTassis fc Mouschovias I ('2005b'), where the magnetic 
shock was shown to be ambipolar-diffusion driven, and 
the associated magnetic cycle was shown to repeat on 
the ambipolar-diffusion timescale. 

4. CONCLUSIONS AND DISCUSSION 

Decoupling of the magnetic field from the matter refers 
to the stage in which the magnetic field becomes dynam- 
ically unimportant for the evolution of the neutral fluid 
and vice versa. We find that the magnetic forces begin 
to decline rapidly with respect to gravity, at a density 
w 10** cm-3 and radius < 10"^i?o ~ lAU. This marks 
the onset of magnetic decoupling. It begins while am- 
bipolar diffusion is still more effective as a flux-reduction 
mechanism than Ohmic dissipation. This is in agree- 
ment with the findings of iDesch fc Mouschovias I (|2001[ ): 
namely, it is ambipolar diffusion that initiates the decou- 
pling of the magnetic field from the matter. 

The field decoupling is not completed until all charged 
species detach from the magnetic field lines. This event 
does not occur in the model cloud since, by the end of 
this typical run, the electrons are still attached to the 
magnetic field lines. However, in Paper III we show that 
there are combinations of free parameters that lead to 
complete decoupling of the electrons at densities « 10*^ 
cm~^. In the present, representative calculation, this 
density is not achieved before the temperature reaches 



At densities nn,c ^ 10^ cm^'^ the dominant charge car- 
riers are the ions and the electrons; they also are the 
dominant electric-current carriers, with the electrons be- 
ing slightly less effective than the ions. In the protostel- 
lar fragment, at densities nn,c > 10* cm~^, the grains 
(negative and positive) are the dominant charge carriers. 
They also become the dominant current carriers above 
?in,c ~ 10*^ cm"'^. 

After the onset of the decoupling of the field from the 
matter, Ohmic dissipation "catches up" with ambipo- 
lar diffusion in effectiveness as a fiux-reduction mech- 
anism. This equality in the ambipolar-diffusion and 
Ohmic-dissipation timescales occurs at a mass-to-flux ra- 
tio w 30 times the critical value. 

The magnetic flux in the innermost part of the model 
cloud, even after Ohmic dissipation becomes dominant, 
is not destroyed. The inward transport of the magnetic 
flux is halted and the magnetic field becomes spatially 
uniform. Indeed, this is exactly the behavior expected 
for a field anchored to the envelope of the cloud and 
supported by external currents, as is the case both in 
nature and in the model cloud. From a mathematical 
point of view, once ambipolar diffusion and Ohmic dis- 
sipation erase spatial variations of the field in the inner- 
most parts of the core, Ohmic dissipation is unable to 
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operate further. This situation may change quaHtatively 
after the magnetic field of the protostar becomes dipolar. 
Such an effect is not observed prior to the recoupling of 
the magnetic field due to thermal ionization; no signif- 
icant pinching of the field lines occurs up to this point 

The magnetic field strength in the innermost 
part of the core at the end of the calcula- 
tion is 0.14 G, comparable to the magnetic field 
of the early solar system as measured in mete- 
orites (iLevv I ll988t __ [Stacev. Lovering. fc ParrTI Il961l : 
iHerndon fc Rowe Ill974( l. 

At the end of the calculation, the hydrostatic core at 
the center of the supercritical protostellar fragment con- 
tains about O.O1M0. The mass, of course, continues to 
increase in time. If one wants to follow the evolution of 
the opaque core until its mass becomes comparable to 
a stellar mass, one would need to include new physics 
in the formulation of the problem (e.g., radiative trans- 
fer, which would allow determination of the temperature 
evolution; thermal ionization; and rotation - see discus- 
sion below) as well as new chemistry. The evolution in 
this new regime will be addressed in a future publication. 
However, the evolution of the isothermal disk surround- 
ing the opaque region until 1 Mq has been accumulated 
at the center has been studied by Tassis & Mouschovias 
(2005a, b). They found that accretion onto the forming 
protostar occurs in a time-dependent, spasmodic fash- 
ion, due to a series of quasi-periodic magnetic shocks. 
In the work discussed here, we have seen the formation 
of the first of this series of shocks, which occurs outside 
the hydrostatic core, as a result of a local concentration 
of magnetic flux.^ Our calculation demonstrated that, 
by the time thermal ionization becomes important, the 

® Note that the creation of the magnetic wall found by Tas- 
sis & Mouschovias (2005b), which is responsible for the series of 
magnetically-driven shocks, is a local effect and does not depend 
on the amount of magnetic ffux accumulated in the central region. 



mass-to- flux ratio of the hydrostatic protostellar core has 
increased by more than two orders of magnitude from 
its critical value. For comparison, observed protostellar 
magnetic fields imply mass-to-flu x ratios two to four o r- 
ders of magnitude above critical (jMouschovias |[T987bf) . 
Although the present result refers to a small protostellar 
mass, the work of Tassis & Mouschovias (2005a, b) has 
demonstrated that ambipolar diffusion alone in the ac- 
cretion phase increases the mass-to-flux ratio of a 1 Mq 
protostar by a similar amount. Interestingly, Shu et al. 
(2006) suggested that in the accretion phase, Ohmic dis- 
sipation alone cannot account for the required flux loss 
despite the fact that, in their semi-analytic approach, the 
magnetic field has no dynamical effect on the system, and 
free-fall accretion is assumed. 

Inside the hydrostatic core, where the thermal-pressure 
force becomes the main opponent of the gravitational 
force, the geometry tends to change from disklike to 
spherical (see, also, Desch & Mouschovias [2001], who 
studied the effect of this change on the results up to den- 
sities of 3 X lO^^cm"'^). The change in geometry changes 
the dependence of the gravitational field from to . 
However, in the mass and density regime we study, this 
effect is quantitatively small. Inside the hydrostatic core, 
the thermal-pressure force reduces the contraction veloc- 
ities to negligible values in any case. At the boundary 
of the hydrostatic core at the end of the run, the differ- 
ence between the gravitational field calculated under the 
disklike and spherical geometry assumptions is only 9%. 
Future calculations pursuing the detailed structure of the 
hydrostatic core at even higher densities will address this 
change in geometry. 

An associated issue is whether the assumption that 
centrifugal forces are negligible is justified. As shown by 
Basu & Mouschovias (1994), by the time a supercritical 
core forms, enough angular momentum has been removed 
by magnetic braking to render the centrifugal force in- 
significant relative to the gravitational force. After this 
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Fig. 12. — (a) Central value of the ratio of ambipolar-diflfusion and Ohmic-dissipation timescales as a function of central density, (b) 
Radial profile of the same ratio at different times, as in Fig. |[2}. 



time, as long as the disklike geometry is preserved, both 
gravity and the centrifugal force scale as r~^, and their 
ratio is therefore constant. However, once the geometry 
changes from disklike to spherical, gravity scales as r"^, 
and a centrifugal barrier can, in principle, be present. 
In our case, the relevant question is whether such a bar- 
rier exists beyond the boundary of the hydrostatic core, 
where the inflow is still finite. Since the deviation of 
the gravitational force from its value calculated assum- 
ing disklike (instead of spherical) geometry is only 9%, 
no significant centrifugal forces are expected. This can 
be quantified easily. From Basu & Mouschovias (1994) 
Fig. 8, we find that the specific angular momentum cor- 
responding to a central mass of 0.01 Mq (the mass of our 
hydrostatic core at the end of the run) in their fiducial 
run is 10^^ cm~^s~^. Under conservation of angular mo- 
mentum, we can calculate the ratio of the centrifugal and 



gravitational forces at the boundary of the hydrostatic 
core assuming spherical geometry. We find this ratio to 
be 6 X 10~^. No realistic cloud rotation can give a sig- 
nificant centrifugal force at the stages of the evolution 
studied in the present investigation. Hence, neglecting 
the centrifugal force is justified throughout each of our 
simulations. 
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